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Abstract Wolf-Rayet stars are the evolved descendents of massive stars. Their 
extraordinary properties make them useful tracers of the stellar initial 
mass function (IMF) in a young stellar population. I discuss how the 
interpretation of spectral diagnostics are complicated by the interplay of 
stellar, nebular, and dust properties. There is mounting observational 
evidence for spatial inhomogeneities in the gas and dust distribution. 
The interplay of these inhomogeneities can significantly alter frequently 
used star-formation and IMF indicators. Specific examples presented in 
this contribution are the starburst galaxies NGC 1614, NGC 2798, and 
NGC 3125. 

1. Background 

Star formation in powerful starbursts is an ubiquitous phenomenon 
both locally and at the highest observable redshifts. A prominent subset 
of the starburst class are the so-called Wolf-Rayet (W-R) galaxies, whose 
observational characteristic is the broad emission bump around 4640 - 
4690 A (Conti 1991). The compilation by Schaerer et al. (1999) lists 
139 members. 

W-R galaxies are important because they permit a study of the star- 
formation properties via stellar spectral features, as opposed to indirect 
tracers based on gas and/or dust emission. Hot stars are notoriously 
elusive even in the strongest starbursts because their spectral signa- 
tures are too weak, coincide with nebular emission lines, or are in the 
satellite-ultraviolet (UV). W-R stars are the only hot, massive stellar 
species detectable at optical wavelengths. This is because they have the 
strongest stellar winds, which in combination with their high tempera- 
tures produce broad (~1000 km s _1 ) emission lines not coinciding with 
emission from H II regions. Examples are N III A4640, C III A4650, 
He II A1640 and A4686, C III A5696, and C IV A5808. The mere de- 
tection of such features proves the presence of stars with masses above 
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40 - 60 Mq (depending on chemical composition) since only stars this 
massive can evolve into the W-R phase. This powerful diagnostic can be 
used for, e.g., inferring a massive star population when the space- UV is 
inaccessible due to dust (infrared [IR]-luminous galaxies), or when broad 
nebular lines veil the O stars (Seyfert2 galaxies). Vice versa, the rela- 
tive strength of the various W-R features gives clues on the distribution 
of W-R subtypes, which in turn provides much sought constraints on 
stellar evolution theory. 

2. W-R stars as proxies for ionizing O stars 

We are currently involved in an optical+near-IR survey of luminous 
starburst galaxies (Leao, Leitherer, & Bresolin, in prep.). Our goals are 
(i) to establish an unbiased sample of metal-rich W-R galaxies, (ii) to 
investigate the stellar content from purely stellar tracers, (iii) to cross- 
calibrate stellar and nebular diagnostics in metal-rich galaxies, (iv) to 
correlate the derived stellar initial mass function with host galaxy pa- 
rameters, and (v) to test stellar evolution theory in high-metallicity 
environments. Our new survey of starburst galaxies is drawn from 
an IRAS sample of Lehnert & Heckman (1995). The galaxies have 
10 < logLiR/L < 11.5, warm IR colors, and are not AGN dominated. 
Detection of the W-R features, together with other standard diagnostics 
allows us to probe the stellar content with a suite of models we have 
developed (Leitherer et al. 1999). In particular, we can search for or 
against evidence of a peculiar IMF at high metallicity, as indicated, e.g., 
by IR observations. 

Thornley et al. (2000) carried out an ISO spectroscopic survey of 27 
starburst galaxies with a range of luminosities from 10 8 to 10 12 L . The 
[Ne III] 15.6 ^m and [Ne II] 12.8 fim lines are particularly useful. The 
ionization potentials of neutral and ionized Ne are 22 and 41 eV, respec- 
tively, the two lines are very close in wavelength (12.8 and 15.6 fJ.ro), and 
they have similar critical densities. This makes the line ratio a sensi- 
tive probe for various star formation parameters, in particular the upper 
mass cutoff of the IMF (M up ) and the age and duration of the starburst. 
Photoionization models of Thornley et al. and Rigby & Rieke (2004) 
suggest that stars more massive than about 35 to 40 M are deficient 
in the observed sample. Either they never formed because of a peculiar 
IMF, or they have already disappeared due to aging effects. This result 
echos that obtained from ground-based near-IR spectroscopy: the strate- 
gic recombination lines are powered by a soft radiation field originating 
from stars less massive than ~40 M (Doyon et al. 1992). An upper- 
mass cutoff as low as 40 M , however, is difficult to reconcile with the 
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ubiquitous evidence of very massive stars with masses of up to 100 M 
in many starburst regions (e.g., Leitherer et al. 1996; Massey & Hunter 
1998; Gonzalez Delgado et al. 2002). Therefore the alternative expla- 
nation, an aged starburst seems more plausible. Under this assumption, 
stars of masses 50 - 100 M Q are initially formed in most galaxies, but 
the starbursts are observed at an epoch when these stars are no longer 
present. This implies that the inferred burst durations must be less than 
a few Myr. Such short burst time scales are surprising, in particular for 
luminous, starburst galaxies whose dynamical time scales can exceed 
tens of Myr. Both the peculiar IMF or the short starburst time scales 
in dusty, IR-bright starbursts are quite unexpected and pose a challenge 
to conventional models in which the starburst is fed by gas inflow to the 
nucleus over tens of Myr as a result of angular momentum loss. 
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Figure 1. Close-up view of the spec- 
tral region around 4650 A in the star- 
burst galaxy NGC 1614. 



Figure 2. Same as Fig. 1, but for 
NGC 2798. Note the broad emission- 
line features from W-R stars. 



Early results from our W-R survey urge caution when relating nebular 
IMF tracers to the actual stellar content. Several metal-rich starburst 
galaxies exhibiting a soft radiation field do in fact have a substantial 
W-R population. As an example, we reproduce in Figs. 1 and 2 portions 
of Keck LRIS spectra around 4650 A of NGC 1614 and NGC 2798, two 
archetypal starburst galaxies with low nebular excitation. We clearly de- 
tect the tell-tale signatures of W-R stars in these galaxies. Unless stellar 
evolution proceeds differently than in our Galaxy, massive O stars must 
be present as well. The fact that we do not directly observe these O stars 
is no contradiction: their spectral lines will be hidden by coinciding neb- 
ular emission lines. We should, however, detect their ionizing radiation 
in a proportion predicted by the measured number of W-R stars and the 
expected ratio of W-R/O stars. The deficit of radiation suggests that 
indirect star- formation tracers, like nebular lines, still require careful 
calibration, in particular when applied to dusty, metal-rich starbursts. 
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After advocating the use of W-R stars to constrain the IMF, we will 
present a cautionary view point in the next section. Interpretation of 
the stellar W-R feature may sometimes be complicated by the inhomo- 
geneous structure of the surrounding interstellar medium (ISM). 

3. The extraordinary W-R cluster in NGC 3125 

As part of a larger project to quantify the stellar and interstellar prop- 
erties of local galaxies undergoing active star formation (Chandar et al., 
in prep.), we have obtained HST Space Telescope Imaging Spectrograph 
(STIS) long-slit far- and near-UV spectra for 15 local starburst galaxies. 
The target galaxies were selected to cover a broad range of morpholo- 
gies, chemical composition, and luminosity. Most importantly, target 
selection was not based on W-R content. 

The UV counterpart of the optical He II A4686 line is the 3^2 
transition of He + at 1640 A. Broad He II A 1640 emission is seen m 
the UV spectra of individual Galactic and Magellanic Cloud W-R stars 
(e.g., Conti & Morris 1990) but is not prevalent in the integrated spec- 
tra of galaxies because of the overwhelming light contribution from OB 
stars in the UV. We observe two different line morphologies of He II 
A1640 emission in local starburst galaxies. Some galaxies show narrow, 
nebular He II emission, while others have a broader profile. If mas- 
sive stars are forming, He II A1640 can sometimes appear as a nebular 
recombination emission line (e.g., Garnett et al. 1991), which has a 
characteristic narrow-line morphology. In the following we will focus on 
one outstanding starburst region in our sample: NGC 3125-1. The mea- 
sured He II FWHM of 5.9 ± 0.4 A is clearly broader than the estimated 
instrumental profile of 2.6 A. The corrected width of the He II A1640 
line is 1000 km s -1 . Similar line widths are typically found in the latest 
subtypes of individual WN stars. 

NGC 3125-1 has by far the largest He II A1640 equivalent width in our 
sample. In Fig. 3 we show the UV spectrum of this object. The 1640 A 
emission is quite prominent and appears significantly stronger than the 
C IV emission at 1550 A. In addition to this feature, we also see N IV 
A1488 and N IV A1720 emission typically found in W-R stars, consistent 
with our interpretation that the strong He II A1640 emission arises in 
the winds of massive stars. We estimate that 6100 WNL stars reside in 
this starburst region, making it the most W-R-rich known example of an 
individual starburst cluster in the local universe. What makes this region 
extraordinary is the number of W-R stars relative to other massive stars. 
The W-R features in Fig. 3 are almost undiluted compared with those 
seen in single W-R stars. Therefore, most of the continuum light must 
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be emitted by the very same W-R stars. Quantitative modeling leads to 
approximately equal W-R and O-star numbers. Such an extreme ratio 
is excluded by stellar evolution models, even for a most unusual IMF. 
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Figure 3. Far-UV spectral region of the starburst region NGC 3125-1. Strong, 
broad He II at 1640 A indicates several thousand W-R stars. The W-R stars provide 
a substantial fraction of the continuum light as well (Chandar et al. 2004). 

How can we reconcile these observations with our current understand- 
ing of W-R stars? The UV continuum slope is normally dominated by 
OB stars; however, in this case the large number of W-R stars also makes 
a significant contribution. One possible explanation is that we are seeing 
the W-R stars through a "hole" where their energetic winds have blown 
out the natal cocoon earlier than have the OB stars. The often large 
reddening derived in starburst regions is consistent with the generally 
accepted scenario in which very young clusters remain embedded in their 
natal material until energetic stellar winds from evolving massive stars 
blow out the surrounding gas and dust. 

If W-R stars are preferentially less attenuated than OB stars in NGC 
3125-1, the equivalent width of He II A1640 and other W-R lines is skewed 
towards artificially large values because for a standard stellar population 
surrounded by a homogeneous ISM, the continuum is emitted by the less 
massive OB stars. If the ISM is inhomogeneous, a textbook property 
of a spectral line associated with a single star changes: the equivalent 
width becomes reddening dependent. 

4. Conclusions 

Determining the massive-star IMF in starburst regions often relies on 
indirect methods, such as the recombination luminosity of the ionized 
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gas. While this technique has been demonstrated to be quite reliable 
in normal galaxies (Kewley et al. 2002), it may break down in extreme 
cases. The conditions prevailing in dusty, metal-rich starbursts may lead 
to a decoupling of the emitted ionizing and the observed recombination 
luminosity. The reasons are still not fully understood but may be re- 
lated to a combination of the destruction of photons by dust and the 
spatial morphology of the dust. The latter effect is particularly rele- 
vant in powerful starbursts where supernovae and stellar winds can have 
large effects on the ISM. Starburst models as well as observations of 
galactic superwinds indicate that 1% to 2% of the bolometric luminosity 
of starbursts is converted into mechanical luminosity. This mechanical 
energy input can remove interstellar gas from the star-formation site on 
a time scale much shorter than the gas consumption time scale, which 
can be hundreds of Myr. The resulting ISM structure invalidates the 
often made assumption of isotropy and homogeneity. As a result, IMF 
indicators can mimic anomalous stellar populations. In some cases, even 
purely stellar indicators may change. If different stellar phases are asso- 
ciated with different dust columns, the stellar equivalent widths in the 
integrated population spectrum may probe dust morphologies, rather 
than the IMF. 
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